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ABSTRACT 

The High Mass X-ray Binary (HMXB) SMC X-l demonstrates an orbital variation of ^3.89 days and a 
super-orbital variation with an average length of ~55 days. As we show here, however, the length of the 
super-orbital cycle varies by almost a factor of two, even across adjacent cycles. To study both the orbital and 
super-orbital variation we utilize lightcurves fro m the Rossi X-ray Tim ing Explorer All Sky Monitor (RXTE- 
ASM). We employ the orbital ephemeris from Wojd owski et al.1 (119981) to obtain the average orbital profile, 
and we show that this profile exhibits complex modulation during non-eclipse phases. Additionally, a very 
interesting "bounceback" in X-ray count rate is seen during mid-orbital eclipse phases, with a softening of the 
emission during these periods. This bounceback has not been previously identified in pointed observations. 
We then define a super-orbital ephemeris (the phase of the super-orbital cycle as a function of date) based on 
the ASM lightcurve and analyze the trend and distribution of super-orbital cycle lengths. SMC X-l exhibits a 
bimodal distribution of these lengths, similar to what has been observed in other systems (e.g., Her X-l), but 
with more dramatic changes in cycle length. There is some hint, but not conclusive evidence, for a dependence 
of the super-orbital cycle length upon the underlying orbital period, as has been observed previously for Her X- 
1 and Cyg X-2. Using our super-orbital ephemeris we are also able to create an average super-orbital profile 
over the 71 observed cycles, for which we witness overall hardening of the spectrum during low count rate 
times. We combine the orbital and super-orbital ephemerides to study the correlation between the orbital and 
super-orbital variations in the system, but find that the ASM lightcurve provides insufficient statistics to draw 
any definitive conclusions on possible correlations. 

Subject headings: accretion, accretion disks - neutron star physics - X-rays:binaries 



1. INTRODUCTION 

■ SMC X-l, firs t discovered with Uhuru observations 
OO ' dLeonget all 1 19711) . is a High Mass X-ray Binary (HMXB) 
£f~) ', cons isting of a neutron star (an X-ray pulsar with 0.7 1 s pe- 
f^*) ■ riod .lLucke et al.ll 9761) and a young B0 supergiant compan- 
O ' ion d Webster et al.lll972t iLillerl 1 1 9731) . An accretion disk is 

, formed around the neutron star, likely partly via wind-fed ac- 

■ cretion where a strong stellar wind from the companion star 
' blows mass beyond its Roche Lobe radius and into the gravi- 

f^i , tational influence of the neutron star. Our view of this system 
J* ■ is at high inclination, as we witness X-ray eclipses of the neu- 

. £h ' tron st ar and disk by the c ompanion star once every orbital 

K> period dSchreier et al.lll972l) . 

These previous studies of the SMC X-l system measured 
' the orbital period to be approximately 3.892 days. As has 
been observed in similar HMXB systems, SMC X-l also 
exhibits a long time scale (sa 60 days) super-orbital varia- 
tion in its X-ray lightcurve dGruber & RothschilJll984l) . Un- 
like the super-orbital va riation in systems such as Her X-l 
dTananbaum et al.lll972[). which has a re latively predictable 
35 day period length (IStaubert et al.lll983l) . the super-orbital 
cycle length in SM C X-l is highly variable and follows no 
obvio us pattern dGruber & Rothschild! 19 84; Woidows ki et al.l 
[T998b . As we elaborate upon further below, lengths of the 
super-orbital cycles in SMC X-l can vary by up to a factor of 
two. 

Accretion disk systems with super-orbital variations have 
been explained with warped disks, seen close to edge-on such 
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that the war p partially obscures our view of the X-ray source 
dKatzlll973l) . Such warps are possibly due to an instability 
driven i n the outer disk by radiation from the central X-ra y 
source (lPettersonlll977t lPringldll996l; iMalonev et al.lll996l) . 
or a numbe r of other mechan i sms (s ee lCaproni et al .120061 for 
a review). Woj dowski et al.1 0998) suggested that the super- 
orbital variation seen in SMC X-l is indeed due to obscuration 
by such a warped disk, as opposed to intrinsic flux variations, 
since the flux and spectrum during orbital eclipse are fairly 
insensitive to whether the system is in a 'low' or 'high' state 
of the super-orbital variation. 

It also has been suggested that SMC X-l possibly has multi- 
ple warp modes in its accretion disk in order to account for the 
wide variation in super-orbital cycle length ( Clarks on et al.l 
2003). For radiatively driven warps, theoretical studies show 
a number of possible "branches" of mode solutions that en- 
compass both prog rade and retrograde, as well as stable 
and d amped warps (Wiiers & Pringle 1999; Qgilv ie & Dubusl 
1200 lb . On the other hand, observational evidence has been 
found in similar systems suggesting interaction between the 
super-orbital variations and the orbital period or som e other 
underlying 'fundamental clock' dBovd & Sm ale 2004). Such 
systems have been shown, in some cases, to have super-orbital 
cycle lengths equal to integer or half-integer multiples of an 
underlying clock, which in the case of Cyg X-2 is the or- 
bital period, while in the case of LMC X-3 and Cyg X-3 is 
not s imply related to any k nown dynamical period in the sys- 
tem (Boyd & Smale 2004). (The' super-orbital variations in 
LMC X-3, howe ver, are likely due to flux variations rather 
than obscuration: IWilms et al.ll200lL ) In Her X-l the super- 
orbital cycle length mainly exhibits three val ues randomly, 
each differing only by half the orbital period (IStaubert et al.l 
[l9l^lStiIi~& Boydll2004tlKlochkov et al.ll2006l) . 

In this work we use data from the All Sky Monitor (ASM) 
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on board the Rossi X-ray Timing Explorer (RXTE) to make 
a comprehensive study of the long term X-ray light curve of 
SMC X-l . We study and characterize both the average orbital 
variation and super-orbital variations in this system. The out- 
line of our paper is as follows. In ^2]we discuss the extraction 
and reduction of the ASM data. We first use these data to 
characterize the average X-ray properties of the orbital period 
(® . Next, in §|4]we discuss how we define the ephemeris for 
the super-orbital variations, and then in §|5]we present their av- 
erage X-ray properties. In §|6] we consider jointly the average 
X-ray properties of the orbital and super-orbital variations. Fi- 
nally, we summarize our conclusions and make comparisons 
to other systems, such as Her X-l, in ^7] 

2. ASM DATA REDUCTION 

The ASM consists of three scanning shadow cameras, each 
equipped with a position sensitive proportional counter that 
views the sky through a set of slits to measure the rela- 
tive intensities and positions of X-ray sources in the sky 
dLevine et al.ll 19961) . The detector operates by comparing ob- 
servations of intensities in any given area of the sky to a cat- 
alog of known X-ray sources to obtain a lightcurve (on time 
scales as short as 90 s) for each target in the field of view. The 
ASM records data in three different energy channels. Chan- 
nels one, two and three are sensitive to X-rays of energy 1 .5- 
3 keV, 3-5 keV and 5-12 keV, respectively. 

We obtained definitive lightcurves for SMC X-l from the 
ASM source catalog at the NASA RXTE website 1 , and an- 
alyzed data beginning MJD 50088.4 and ending 54083.8. 
As mentioned previously, the ASM data consists of a least- 
squares fit to the sky's modeled X-ray spectrum for all the 
known sources within its field of view. For this reason the 
count rate (cps) measurements from the satellite will occa- 
sionally fall below zero when the target source counts are 
lower than the difference between the expected background 
and the actual background; we have not excluded such nega- 
tive counts from our extracted lightcurves. 

Our analysis was perfo rmed using ISIS version 1.4.2-5 
dHouck & Denicolal [2000), using custom scripts written in 
S-lang 2 (the scripting language embedded within ISIS), as 
well as routines publicly available from the S-lang/ISIS 
Timing Analysis Routines 3 (SITAR). Using these 
routines to read the ASM data, we only retained data points 
where the ASM solution had a maximum x 2 -value of < 1.5. 
Before beginning any analysis, we performed a barycenter 
correction on the lightcurve to account for the movement of 
the satellite within the Earth-Sun system. 

The ISIS def ine_counts function allows arbitrary ar- 
rays of lo/hi bin, value per bin, and error on that value to be 
registered as a finable data set, with presumed diagonal re- 
sponse matrix and unity effective area. We used this function- 
ality to register the folded ASM lightcurves of $3] and SI for 
fitting, and used combinations 4 of custom defined S-lang 
functions to fit these data. When it was necessary to fit a 
function to (x,y) pairs rather than histogram data sets (see 

1 http://heasarc.nasa. gov/ doc s/xte/ ASM/ sources . html 
2 www . s— lang . org 

3 |http: // space .mit ■ edu/CXC/analysis/SITAR| 
4 As opposed to IDL, MATLAB, or ISIS's own array-based fitting - where 
each unique model must be a 'separate' function - ISIS histogram fits use 
S-lang to parse model expressions such that almost any 'mathematically 
sensible' combination of individual model components results in a properly 
defined total model. (This functionality will soon be extended to the ISIS 
array-based fitting as well; J. Houck, priv. comm.) 
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FIG. 1. — The ASM lightcurve of SMC X-l folded on the orbital period, 
fit with a 10 parameter empirical model. Error bars are 1 cr, and are derived 
from the standard deviation of the mean rate in each bin. The bounceback 
feature is clearly visible between phases .875 and 1.025 

we used the ISIS function f it_array, which takes as in- 
put arrays of x- and y-values, the relative weights of the data 
points, parameter initial values and limits, and a reference to 
a single user-defined S-lang fit function. 

3. FOLDING THE ORBITAL PERIOD 

Our first step in characterizing the behavior of SMC X- 
1 was to search the ASM lightcurve for evidence of the 
aforement ioned periodicities. Using a Lomb-Scargle peri- 
odogram (iLombll 19761: IScargldll982l ) we were able to detect 
both a 55 day period (the average length of the super-orbital 
cycles; see © as well as a period of approximately 3.89 
days (i.e., the orbital period). These findings correspond to 
previous observ ations of the source (Wojdowskiet alj|1998t 
I Wen et ai1l2006l). We obtaine d an o rbital ephemeris for the 
source from IWoidowski et al.l dl998l) . which was defined in 
terms of the initial observation time, the period and the period 
first derivative. Using this ephemeris with the S I TAR routine 
sitar_pf old_rate, we folded the ASM lightcurve into 
60 phase bins, as shown in Fig.[T] 

We had originally thought to search for improvements on 
the Wojdowski et al. ephemeris by detecting a drift in the 
zero point of orbital phase (even relative to the known or- 
bital period derivative from that work) over the 1 1 years of 
ASM data we analyzed, but found that any possible drift in 
the ephem eris was well with i n the observational errors of the 
ASM. The Woidows ki et al.l (|T998) ephemeris remains accu- 
rate. 

Using the folded lightcurve, we were able to observe sev- 
eral features of the orbital profile. The most interesting of 
these features is a "bounceback" phenomenon occurring dur- 
ing the eclipse. Because we are observing the system at high 
inclination, we expect to have few to no counts during the 
eclipse phases, when the companion star obscures the line of 
sight between the satellite and the neutron star-disk system. 
Surprisingly, we observe the low counts expected on either 
side of the eclipse, but witness a small rise and fall of counts 
during the mid-eclipse phases. We will henceforth refer to this 
region of recovering count rate as the "bounceback." Another 
interesting feature is a complex modulation, modeled below 
as sinusoidal, that occurs throughout the signal profile, and 
that is very visible during the non-eclipse phases. 

We were able to obtain a fit to describe these features of the 
orbital profile using a ten-parameter empirical function. Our 
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FIG. 2. — Folded orbital profiles of the three ASM energy channels, from left to right 1.5-3, 3-5 and 5-12keV. Both the bounceback and complex orbital 
modulation are more prominent in the lower energy channels. 



final model was obtained as follows. We modeled the bounce- 
back feature in the center of the eclipse with a parabola. Al- 
lowing a free parameter for the magnitude of the curvature, 
which could range to any real value (positive or negative), 
we noted that the curvature was fit to positive values in all 
channels. Eclipse ingress and egress were modeled by lin- 
ear functions, dropping from/rising to a constant value out- 
side of eclipse. We multiplied the entire function by a sine 
wave plus a constant to account for the modulation observed 
in the profile, and subtracted a Gaussian from the center to ac- 
count for a broader component of the eclipse, possibly due to 
partial obscuration of potentially extended X-ray emission by 
the companion star and/or accretion stream. The center of the 
Gaussian and parabola were fixed at 1, while the sinewave fre- 
quency was fixed at an integer or half-integer multiple of the 
inverse orbital period. We attempted the fit with several dif- 
ferent integer and half-integer values and used the frequency 
with the lowest \ 2 m our final fit function. It is interesting to 
note that the best fit was obtained modeling the modulation 
with a sine wave of frequency 4 (in units of the inverse orbital 
period), implying a complex disk warp structure. 

Aside from the amplitude of the constant, amplitude and 
phase of the sine wave, and area and standard deviation of 
the Gaussian, free parameters also included widths for each 
of the regions: the eclipse, the ingress, and the egress. The 
folded data, fit and residuals are plotted in Fig. Q] and the pa- 
rameters are given in Table 1. This model, as well as other 
non-standard models used for fitting in this paper, were writ- 
ten by the authors and can be made available upon request. 

Once a successful model was obtained we were able to uti- 
lize it to fit orbital profiles to each of the individual energy 
channels of the ASM data (1.5-3, 3-5, and 5-12keV) and 
compare the normalizations of features of the system in dif- 
ferent energy ranges. Fig. [2] shows profiles and fits in each 
of these three energy channels, and the fit parameters are also 
presented in Table 1 . In general we found that the lightcurve 
count rate was dominated by the high energy portion of the 
spectrum. As expected we also found more modulation in the 
low energy channels than in the high, likely due to preferential 
absorption of low energy X-rays, and also possibly due to the 
relative sizes of the sources of high and low energy X-rays. 

Due to count rate uncertainties, it is not initially clear 
whether or not the bounceback and sinusoidal variation are 
significant, real features of the orbital profile. In order to fur- 
ther explore the significance of the bounceback, we performed 



fits of the orbital profile without the parabolic term, resulting 
in a Ax 2 from the fits to the full model described above of 
14.1 in the summed channel fit and 1.0, 2.0 and 10.0 in fits 
of the first, second and third channels respectively. We also 
performed the orbital profile fits without the sinusoidal mod- 
ulation, which resulted in a A\ 2 of 13.2 in the summed chan- 
nel fit and 5.4, 7.9 and 2.3 in the fits of the first second and 
third channels respectively from the x 2 -values obtained using 
the full model. The statistics are not good enough to allow 
study of the bounceback or the sinusoidal modulation in in- 
dividual periods. The bounceback is not clearly exhibited in 
averaging of most short portions of the lightcurve (i.e., only 
a few orbital periods) but can only be found with the statis- 
tics provided by long averaging times. A more thorough un- 
derstanding of the bounceback would require better statistics 
than presently available from the ASM data. 

In order to get a more quantitative idea of the energy distri- 
bution of the emission, we created a folded color lightcurve, 
where we defined the color, C, as 

cM-ch3 

chl+ch3 

Here chl is the count rate measured in ASM channel one, and 
ch3 is the count rate measured in channel three. Thus if the 
emission at a certain phase is dominated by hard X-rays, then 
the value of the color will be negative, and if the emission 
is dominated by soft X-rays the color will be positive. The 
orbital phase-dependent color lightcurve (with 40 phase bins) 
is plotted in Fig. [3] 

Overall the color lightcurve demonstrates that the source is 
relatively monochromatic and dominated, as observed in the 
channel fits, by hard X-ray emission. We did find that at the 
edges of the eclipse, where the source was being partially ob- 
scured by the companion star, that the X-ray emission possi- 
bly becomes much harder. This might be attributable to soft 
X-ray absorption from the outer atmosphere of the compan- 
ion star; however, this region of the lightcurve is also where 
the average chl lightcurve dips below 0. If we apply a sys- 
tematic shift upward of 0.05 cps to each of the lightcurves, 
this hardening disappears, although the hardening remains a 
possibility within the uncertainties. On the other hand, the 
bounceback emission, in contrast to the emission in the rest 
of the orbital cycle, clearly becomes softer. 

4. DEFINING THE SUPER-ORBITAL EPHEMERIS 

The inherent noise of the lightcurve and the highly vari- 
able nature of the super-orbital cycle length provided some- 
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TABLE 1 

Orbital Profile: Constant*((1.0+Sinewave)*Orbit+Gauss) 



ASM 




Constant 


Sinewave 






Orbit 






Gauss 


Channel 


(50 DoF) 




Norm 


Phase 


Lwidth 


Rwidth 


Cwidth 


Height 


PNorm 


Area 


Sigma 


1+2+3 
1 
2 

3 


1.00 
0.95 
0.78 
1.14 


1 74+O.O2 
40+° 05 
0.39 + $ 
0-86^ 


036+ 016 
UUJO -0.014 

0.063!^f 
0.048^;0|6 

o.oi8i8:8!^ 


0.68;»^ 
n 70+0.09 

U/S -0.09 
74+0-26 
u - '^-0.74 


n mo 40 001 

uluu -0.005 
i 17 +0.017 
u,11 '-0.015 
(1 1 00+0006 

0QS+0 002 
u - uy;, -0.004 


, n 3+0.005 
u - luJ -0.003 

1 12+0-013 

u ' llz -0.019 
106+0 007 

ioo+o .006 


1 34+fl.OOl 
u - 1JH -o.ooi 

122+0- 878 

-0.025 
140+0-860 
u - lw -0.009 

u iJ -0.005 


77 +0.02 
u -"-o.oi 

0.73*}* 

o-87!8:!J 


io+ 03 

u - lu -0.03 

17+0-17 
u.i ;_ Q 47 

12+0- 14 
u - lz -0.09 

09+0 .07 


-4.9!i; 3 8 
-6 3 +1 2 

—2.1 


U1 °-0.02 
41+0.25 

17+0-17 
ul/ -0.06 

ul '-0.04 



NOTE. — ORBIT function parameters: Lwidth=width of the eclipse ingress, Cwidth=widfh of the eclipse from low point on the left of the bounceback to low point on 
the right of the bounceback, Rwidfh=width of eclipse egress. Height=Distance from lowest point in the eclipse to the constant in non-eclipse phases. PNorm=Height of the 
parabola lit in the bounceback (center of the eclipse) region. In Lwidth and Rwidth regions the profile was modeled with a linear function and in the Cwidth region the profile 
was modelled as a parabola. Center of Orbit model and Gauss were frozen at phase=l and the sinewave frequency was frozen at f=4 (Orbital Period) -1 . Error bars are 90% 
confidence levels. 




Orbital Phase 

FIG. 3. — Top: The orbital-folded color lightcurve, exhibiting overall hard 
X-ray dominated emission with softer X-ray emission in the bounceback re- 
gion. Error bars are are 1 a error bars. Bottom: The same folded lightcurve 
as above; however, each energy channel lightcurve was systematically shifted 
upward by 0.05 cps. 

what of a challenge as far as the definition of the super- 
orbital ephemeris (super- orbital phase to date correlation) 
was concerned. Here we wish to consider only variations on 
the longer super-orbital time scale, and not those produced 
by the orbital variations. We therefore removed points that 
fell within the eclipse of the neutron star by the compan- 
ion from the lightcurve. Using our orbital fold, discussed 
in the previous section, we determined the orbital phases of 
the eclipse to lie between 0.85 and 1.15 (Fig. [TJ, and used 
the IWoidowski et alj (1 19981) ephemeris to calculate the time 
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FIG. 4.— SMCX-1 lightcurve between MJD 50150 and 50475, smoothed 
with a Gaussian filter. Fourth order polynomial fits in the minimum regions 
are superimposed. Minima of fits taken as the zero phase points for the super 
orbital period are denoted by vertical lines. 

ranges during which these phases would occur. We then re- 
moved all ASM data points that fell within those ranges. 

Using the lightcurve with the eclipse phases removed, the 
ASM lightcurve remains intrinsically noisy. For example, the 
peaks of each super-orbital cycle have a number of promi- 
nent dips (see Fig. |4}, which appear to be akin to the "pre- 
ecl ipse dips" seen in the l ong term X-ray ligh t curves of Her X- 
1 (IShakuraet all fl998l: IStelzer et all Il999l: iKlochkov et all 
2006). To further reduce the inherent noise in the ASM 
lightcurves of SMC X-l, we applied a simple Gaussian 
Smoothing to the unbinned ASM data via an FFT 5 with a 
smoothing scale of 16 data bins (~ 2.5 days on average). 
We defined low count regions by selecting regions of the 
lightcurve that were below a cutoff count value (< 1.675 cps), 
and that were wide enough (> 7 days) to assure that they were 
not residual noise fluctuations from the lightcurve. We then 
fit fourth order polynomials to each of the low count regions 
using the ISIS function array_f it, choosing a uniform 
weighting of the data points (see We completed the def- 
inition of the super-orbital ephemeris by defining the minima 
of the polynomials in the low count regions as points of zero 
phase. This procedure is illustrated in Fig. |4] which shows a 
portion of the smoothed lightcurve with the polynomials and 

5 Although only strictly valid for an evenly spaced lightcurve, this proce- 
dure created a well-behaved lightcurve, comparable to various schemes with 
wweighted averages, or binned schemes with padding for data gaps, that we 
also tried. 
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FIG. 5. — The time evolution of super-orbital cycle length. Error bars repre- 
sent the uncertainty in cycle length due to the uncertainty in the determination 
of super-orbital zero phase times. The oscillatory structure superimposed on 
the long- term trend has previo usly been attributed to multiple warp modes in 
the disk jClarkson et alJ2003l) . 
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FIG. 6. — Dynamical Lomb-Scargle periodogram for the entire ASM 
lightcurve of SMC X-l. The periodogram is calculated in 315 day intervals, 
slid 50 days each bin. (See lWilms et aU200U ') 

minima (zero-phase points) superimposed. 

With the ephemeris defined as in Fig. @] we were able to 
extract the start and stop time of each super-orbital cycle. We 
estimate, based upon multiple trials with variants of our zero- 
phase search routines (weighted averages, binned routines, to 
define the lightcurve, and several different functional forms 
to fit the minima) that our cycle start times have an accuracy 
of ±2 days. To search for possible patterns, we plotted the 
length of each super-orbital cycle versus the average date dur- 
ing that cycle. This plot is shown in Fig. [5] There appears 
to be a somewhat oscillatory trend, superimposed on a much 
longer time scale pattern. The longer term trends, especially 
the short (40-50 day) super-orbital cycles near the beginning 
of the ASM lightcurve, are likely what have been identifie d 
as multiple warp modes in the disk (i.e., Clar kson et alj 2003). 
Specifically, if we apply a sliding Lomb-Scargle periodogram 
to the ASM data (exactly as w e have previously shown with 
ASM lightcurves of LMC X-3: IWilms eFallllO Ol). the rather 
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FIG. 7. — Histogram distribution of super-orbital cycle lengths with bi- 
modal (sum of two Gaussians) fit superi mpose d. The fit shown in this figure 
was used in the simulations described in 36.21 

ragged pattern of Fig. [5]appears more as a smooth evolution of 
super-orbital period, as shown in Fig. [6] Comparing to Fig. [5] 
however, we see that the smoothness of this trend is partly 
an artifact of the averaging process, and that a sliding Lomb- 
Scargle periodogram by itself is unrevealing for the fine detail 
of the evolution of this source. 

We studied the distribution of these varying cycle lengths 
by creating a histogram of their values. The number of bins 
in the histogram was chosen because it produced the most 
easily seen pattern, but does not significantly change the re- 
sults of the analysis. For example, doubling the number of 
bins does not change the shape of the distribution, except to 
add more noise to the overall trend. We found in this case 
that the distribution has a distinct double-peaked pattern, sim- 
ilar to other systems of this type. Specifically, such systems 
as Cyg X-2, LMC X-3, Cyg X-3 and Her X-l show dou- 
bly peaked distributi ons for histograms o f wait times between 
successive minima (Bovd & Smale 2004). It is possible that 
the SMC X-l system is similarly oscillating (somewhat ran- 
domly) between two extreme periods. The profile of super- 
orbital cycle lengths is shown in Fig. [7] 

Once we had defined the ephemeris, we attempted to find 
some correlation between the length of the super-orbital cy- 
cles and the X-ray flux from the system. Towards this end, 
for each super-orbital cycle we binned the rates into a his- 
togram and attempted fits with several different functional 
forms; none of these were successful for all of the super- 
orbital cycles. This was due both to the large variance in 
the shapes of each super-orbital variation (including the afore- 
mentioned dips), as well as the lack of data for some portions 
of the lightcurve. Instead, we performed five point spline fits 
on each of the super-orbital cycles. From the spline fits we ex- 
tracted the minimum, peak and mean flux during each super- 
orbital cycle and plotted each of those values against both the 
super-orbital cycle length, and the length minus its mean value 
of approximately 54.4 days. We were not able to find any cor- 
relations between super-orbital cycle length and flux with this 
method, nor by working on a coarser time scale, i.e., by aver- 
aging over five consecutive super-orbital cycles and searching 
for correlations between the same variables. 

5. FOLDING ON THE SUPER-ORBITAL PERIOD 
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TABLE 2 

Super-Orbital Zero Phase Times and Lengths 



Cycle Number 


Start Time 


Length 


Cycle Number 


Start Time 


Length 


Cycle Number 


Start Time 


Length 




(MJD) 


(Days) 




(MJD) 


(Days) 




(MJD) 


(Days) 


j 


sni 9t s 

JU17J.O 


61 6 


A J 


5 1420 4 


59 6 


4Q 


52824 5 


48 6 


2 


50955 4 


58 2 


26 


51479 9 


69 7 


50 


52873 2 


48 6 




5031 3 6 


56.5 


27 


51549 6 


61 


5 1 


52921 8 


47.8 


4 


50370 1 


50 


28 


51610 6 


63 


52 


52969 6 


61 4 


5 


50420 1 


53.2 


29 


51673 6 


70 1 


53 


53031 


52 5 


g 


50473 4 


46.8 


30 


51743 7 


57 1 


54 


530X3 5 


52 7 


7 


50590 7 


52 


^ j 


5 1 son 8 


56 8 


55 


53136 2 


54 




50572 2 


45 2 


32 


5 1 857 5 

JIOJ/.J 


48 8 


56 


53190 1 


62 3 


9 


50617 4 


47 3 


33 


51906 3 


58 6 


57 


53252 4 


59 6 


10 


50664 7 


40.3 


34 


51964 9 


54 8 


58 


53312 


61.2 


1 1 


50705 


43.0 


35 


52019 7 


66.7 


59 


53373 2 


46.5 


12 


50747 9 


49.9 


36 


590R6 4 


48 8 


60 


53419 7 


66.6 


13 


50797.9 


39.6 


37 


52135.2 


59.7 


61 


53486.3 


62.7 


14 


50837.5 


50.7 


38 


52194.9 


50.6 


62 


53548.9 


59.6 


15 


50888.2 


44.3 


39 


52245.4 


64.8 


63 


53608.5 


62.2 


16 


50932.4 


49.4 


40 


52310.3 


62.3 


64 


53670.7 


54.8 


17 


50981.9 


45.2 


41 


52372.5 


46.8 


65 


53725.5 


54.6 


18 


51027.0 


65.7 


42 


52419.4 


56.6 


66 


53780.1 


54.4 


19 


51092.7 


52.6 


43 


52476.0 


72.1 


67 


53834.5 


44.4 


20 


51145.3 


57.9 


44 


52548.1 


50.5 


68 


53878.9 


47.1 


21 


51203.2 


58.5 


45 


52598.5 


49.1 


69 


53926.0 


47.8 


22 


51261.8 


52.3 


46 


52647.6 


44.1 


70 


53973.8 


42.2 


23 


51314.0 


58.5 


47 


52691.7 


71.9 


71 


54016.0 


39.0 


24 


51372.5 


47.8 


48 


52763.6 


60.9 









NOTE. — We estimate our methods for determining zero phase of the super-orbital cycles yields an uncertainty of about 2 days, which gives 
a corresponding uncertainty of 2 \fl days in the determination of the super-orbital cycle lengths. 





Super-Orbital Phase 

FIG. 8. — Lightcurve Folded on super orbital period with fit of two super- 
imposed Weibull functions. 

With the super-orbital ephemeris defined, we were able to 
assign a super-orbital phase to each point on the lightcurve 
and proceed to fold the data over all super-orbital cycles. 
We chose to average the ASM lightcurve in 25 phase bins 
of the super-orbital cycles, and show the results in Fig. [8] 
It is interesting to note the asymmetry between the rise and 
fall of the curve, with the leading edge (super-orbital eclipse 
egress) exhibiting a much sharper rise than the falling edge's 
decline. The best fit to this profile was obtained with the 
sum of two scaled Weibull functions ( i.e., oc x 1 exp[-(x/ /?)"]; 
see Appendix B of Nowa ket alJ 19991 for historical references 
and description of the uses of this distribution function). We 
elected to use Weibull functions because their asymmetry can 
reasonably describe the asymmetry of the folded lightcurve. 
We used the sum of two curves as this matched the broadly 
peaked nature of the folded profile much better than did a sin- 



0.4 0.6 

Super-Orbital Phase 

FIG. 9. — Color diagram for super-orbital profile. Emission is relatively 
monochromatic and dominated by hard X-rays. In lower intensity regions the 
emission is harder than in high intensity regions, which is in contrast to the 
eclipse emission in the average orbital profile. 

gle function. This was also suggested by the shape of the 
individual super-orbital cycle profiles, many of which were 
dramatically double peaked. Fit results are presented in Table 
3. 

As in our analysis of the orbital profile, using the data from 
the individual ASM energy channels we were able to create a 
super-orbital profile for each of the three energy channels. We 
fit each of these profiles with the same functional form used 
to fit the total profile. Fit results for the individual channels 
are also presented in Table 3. As expected from our analysis 
of the orbital profiles we once again found that the profile was 
dominated by high energy emission. Similar to our analysis 
of the average energy profile of the orbital variation, we cre- 
ated a color lightcurve for the super-orbital profile by plotting 
the color versus super-orbital phase, again using the definition 
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TABLE 3 

Two Weibull Function Fits to Lightcurves Folded on Super-Orbital Cycle 



ASM 




Normi 


Peak, 




ft 


Norm2 


Peak 2 


a-7. 
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(17 DoF) 


















1+2+3 
1 

2 
3 


1.95 
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I2 ml3 
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i 97+0.01 

^'-o.oi 


9 34+0.21 
J -O.02 

2 1 +33 

z,1 -0.8 
2 5 +0 - 7 

2 03+ 015 


26 +a02 
lfi 40 - 17 

U - 1D -0.07 
n 90+O.O6 

20+ 03 
u - -0.01 


1 30+ 009 

1 ' ;>u -0.10 

n 70+0.03 
27+ 02 

u - z '-0.03 

69 +0 00 

u -°'-0.40 


1 57+ 001 

J -0.00 
i ^9+0.04 
1 - JZ -0.03 

1 5Q+0.01 

Jy -0.00 
i C4+O.O2 

'■ JH -o.oi 


4 76+ 09 

o 9C+I.8I 
J - ZJ -0.03 

6 1+ 3 ' 9 

D.l_j 9 


n ss+o.09 
63 +0 33 
0.99!°; 3 > 



NOTE. — Weibull functions had the form Normi(A -A'o)exp[-((x-jco)//3i) ai ], etc., where *o was determined by fixing the location of 
the function maximum to Peak] . Error bars are 90% confidence levels. 



of eq. (1). We present these results in Fig. [9] Once again we 
can see that the emission is fairly monochromatic, with two 
exceptions. Near phase 0.9 of the super-orbital period, there 
is weak evidence for a softening. More clear, however, is that 
near phase of the super-orbital cycle (i.e., the low count rate 
region in which the neutron star is likely eclipsed by the disk) 
the emission appears to become harder, although the color val- 
ues in this region have relatively large uncertainty. This is in 
marked contrast to the "bounceback" region of the orbital pe- 
riod fold, where the emission became softer. As in our color 
plot for the orbital profile, we have examined the same plot 
with a systematic shift upwards of 0.05 cps in all channels, 
however, in this case it does not significantly change the re- 
sult. 

6. JOINT ORBITAL AND SUPER-ORBITAL PROFILES 

6.1. Two Dimensional Phase Folds 

To better understand the link between the orbital and super- 
orbital variations in SMC X-l, we created a two dimensional 
fold of the SMC X-l lightcurve. Once we had defined an 
ephemeris for both the orbital and the super-orbital variations 
(as discussed in $3] and §|4] respectively) we were able to as- 
sign every point on the lightcurve both an orbital and a super- 
orbital phase. We then created a grid of histogram bins with 
orbital phase on one axis and super-orbital phase on the other, 
and sorted the data points into that grid, averaging the inten- 
sity measurements in each bin. The result was the two di- 
mensional histogram shown in Fig.[l0] While we can observe 
the asymmetry noted in the super-orbital profile in this his- 
togram, it does not contain sufficient statistics to give specific 
insight as to the evolution of orbital profile features at differ- 
ent super-orbital phases and vice versa. The choice of binning 
is made for maximum clarity. Unfortunately, current statistics 
are not sufficient such that finer binning would provide more 
information. 

To gain another visualization of the form of the orbital pro- 
file at different super-orbital phases, we also created several 
orbital profiles, each containing only points of the lightcurve 
in a small range of super-orbital phase, and plotted these pro- 
files on a single set of axes. We present these orbital phase 
histograms in Fig. QT] Here again we find that we have 
insufficient statistics to demonstrate any definitive trends in 
the shape of the orbital profile dependent upon super-orbital 
phase. Further study of this system with better statistics (e.g., 
with a series of pointed X-ray observations) are likely neces- 
sary to witness any trends. 

6.2. Searching for Dependencies upon Orbital Period 

In accordance with patterns seen in other sources, e.g., 
Her X-l where the 'turn on' of the super-orbital cycle is tied 




12 3 

Average ASM Rate (cps) 




FIG. 10. — Two dimensional folded lightcurve. Variation in orbital phase is 
plotted along the jc-axis while super-orbital phase is plotted along the y-axis. 
The top panel shows the color scale. 

to a half integer multiple of the orbital period (Staube rt et al.l 
1 1983b iKlochkov et al.ll2006l) or Cyg X-2 where the time be- 
tween s uccessive minima is a n integer multiple of the orbital 
period (iBovd & Smale 2004), we searched for a correlation 
between the orbital period length and the super-orbital cycle 
lengths. We used a set of ^3500 values between one quar- 
ter of the orbital period and one half of the longest observed 
super-orbital cycle to search for a value that divided most 
evenly into all of the super-orbital cycle lengths. For any 
given member of the set, we tested its relation to the super- 
orbital cycle lengths by dividing each of the 7 1 observed cycle 
lengths by that value, and summing the non-integer remain- 
ders of each of the 71 quotients. A value that fit relatively 
evenly into the set of super-orbital cycle lengths was one that 
exhibited a minimum value of the sum of non-integer parts. 
We found 16 values that exhibited minima in this formula, 
which we define as a point for which the sum of non-integer 
parts was <31. (The maximum value of the sum of the re- 
mainders was ~43.6, while the mean of the sum was ~36.2.) 
Of the 16 minima we found 7 that had some integer or quarter- 
integer relation to the orbital period, i.e., approximately 0.25, 
0.5, 1, 2, 3, 5 and 6 times the orbital period length. The other 
9 values that exhibited minima had no obvious relation to the 
orbital period length. 
To test the significance of these observed values, we created 
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Orbital Phase 

FIG. 1 1 . — Lightcurve folded on orbital period for specific ranges of su- 
per orbital phase, demonstrating the correlation between the two. Statistics 
are not good enough to enable identification and analysis of the bounceback 
feature in these folds. 

a histogram of the sum of non-integer parts in order to study 
its distribution. The small value of this sum that we had ob- 
tained for 0.25, 0.5, 1, 2, 3, 5 and 6 times the orbital period 
length appeared rare, but not necessarily unusually so; there- 
fore, we performed a Monte Carlo simulation to determine 
their significance. We first fit a bimodal (sum of two Gaus- 
sians) distribution to the super-orbital cycle length data shown 
in Fig. [7] from which we could randomly choose super-orbital 
period lengths. A single simulation consisted of randomly 
choosing 7 1 super-orbital cycle lengths from this distribution, 
and then performing the same analysis as done for the actual 
data. Specifically, the main goal was to create a histogram 
of the sum of non-integer parts from the quotients of each of 
the 3500 trial values with the chosen set of super-orbital cycle 
lengths. We performed 10 6 such simulations and averaged all 
of the resulting histograms to obtain a theoretical distribution. 

From this simulated distribution we could estimate the like- 
lihood of obtaining the sums of non-integer remainders with 
values less than or equal to those that had been found for the 
quarter-integer and integer multiples of the orbital period. We 
found these probabilities to be *4.2%, 0.1%, 4.2%, 3.6%, 
0.6%, 4.2% and 2.5% for the members of the set that were 
0.25, 0.5, 1, 2, 3, 5 and 6 times the orbital period length, re- 
spectively. From these simulations we concluded that while 
these values appear to have interesting relations to the orbital 
period, and their corresponding sums of non-integer remain- 
ders are somewhat rare, they are not rare enough to point con- 
clusively towards a relationship to the orbital period. Further 
study of this correlation - either via ASM observations with 
a longer baseline, or via pointed observations that can more 
accurately determine the start and stop times of each super- 
orbital period - may prove useful in determining if there truly 
is a relationship between the super-orbital cycle and orbital 



period lengths. 

7. SUMMARY 

In our analysis of the ASM lightcurve of SMC X-l we have 
observed several interesting features of both the orbital and 
super-orbital pro f iles. By using the ephemeris provided in 
Woi dowski etail ([1998), we were able to fold ~ 1 1 years of 
data on orbital period to obtain an average orbital profile. The 
main results of this analysis are as follows. 

• The lightcurve is characterized by complex modulation 
during the persistent emission (here modeled as a sinu- 
soid, but which could possibly fit another model within 
the uncertainty), with a broad, asymmetric eclipse. 

• We were able to witness a count rate bounceback phe- 
nomenon when averaging over large numbers of orbital 
cycles, which consists of a small rise and fall in counts 
during orbital mid-eclipse phases. 

• The emission becomes somewhat softer during the 
eclipse/bounceback, while the persistent emission is 
fairly monochromatic. 

Despite several pointed observations during orbital eclipse, 
in both low and high states of the super-o rbital period, the 
bounceback has not b een previously seen (IWoo etal]ll995t 
Woidowski et al.|[l998l) . We did not have sufficient statistics 
to observe the dependence of this phenomenon on the super- 
orbital cycle, it is therefore possible that the bounceback is 
only present in certain super-orbital phases, or that its pres- 
ence is only intermittent. Given the relatively few pointed 
observations of SMC X-l during eclipse, it is possible prior 
pointed observations have simply not observed the state in 
which the bounceback is most prominent. 

The softening during the eclipse might be akin to that 
seen in so-called A ccretion Disk Corona (ADC) sources 
dWhite & Holt! |1982|) . Such sources show broad, partial 
eclipses, with deeper eclipses in the hard X-ray (for example , 
see the lightcurves of X1822-371. iHeinz & Nowa3l2001h . 
One model of these sources has a very spatially extended 
corona scattering X-ray from the central source into our line 
of sight; however, X-rays traveling through the atmosphere 
just above the midpl ane of the disk are pr eferentially absorbed 
in the soft X-rays dWhite & Holtlll982l) . Thus, the eclipse 
of the base of this extended corona is preferentially blocking 
harder X-rays, and leads to an overall softening of the spec- 
trum. However, if such is occurring in SMC X-l, this does 
not explain the "bounceback" of the lightcurve. Perhaps addi- 
tional scattering into our line of sight at orbital phase 0, from 
the atmosphere of the companion, rather than just an ADC, 
must also be invoked to explain the bounceback. 

The super-orbital variation in SMC X-l is not as easily dealt 
with as the orbital variation because of the complicated nature 
of the super-orbital ephemeris. The super-orbital cycle length 
is highly variable and Lomb-Sc argle techniques alone are not 
sufficient to define its structure (I Wen et al.l2006l) . In our anal- 
ysis we were able to obtain a super-orbital ephemeris from the 
ASM data by modeling the locations of zero phase for each 
of 71 cycles. From that ephemeris we were able to extract 
an average super-orbital profile, as well as information about 
trends in the super-orbital cycle length. Our main results are 
as follows: 
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• In general we observe an asymmetry in the average 
super-orbital profile, most notable in a sharper rise in 
the low state egress than fall in the low state ingress. 

• In contrast to the trend in the orbital profile, where we 
see softer emission during low count states, we find 
that emission appears to harden during low states in the 
super-orbital profile. 

• We have found the lengths of the super-orbital cycles 
in S MC X-l to be b i modal ly distributed, in agreement 
with |Bovd & Smale d2004f) . and have also found sug- 
gestive, although not definitive, evidence that the varia- 
tion in length is driven by some internal clock, the fun- 
damental period of which is related to th e orbital period 
as ob served in Her X-l and Cyg X-2. (Bovd & Smale 
l200l . 

• In contrast to fin dings in analysis of Her X-l 
(IStill & BoydH 2004). we have been unable to find any 
correlation between the super-orbital cycle length, and 
X-ray flux. 

• In contrast to IClarkson et alJ ((2003), who report a 
smooth evolution of the super-orbital cycle length, we 
observe sharp variations, often with large changes in 
length between adjacent super-orbital cycles. 

This latter result is especially evident in observed super- 
orbital cycles 41 through 44 (MJD 52372.5 to 52598.5), in 
which the cycle duration rises from —47 days in cycle 41 to 
—72 days in cycle 43, and falls back to —50 days in cycle 44. 
Likewise between super-orbital cycles 46 through 49 (MJD 
52647.6 to 52873.2) there is a rise from -44 days to -72 
day s in cycles 46 and 47, and a fall back to —49 days in cycle 
49. IClarkson et al.l (|2003) also theorizes that the evolution of 
super-orbital cycle length is possibly sinusoidal. The lack of 
a regular, sinusoidal pattern in our data can easily be seen in 
Figs.|5]and|6] 

In our spectral analysis of 7 1 super-orbital cycles we were 
able to observe a clear super-orbital profile in the first channel 
(1.2-3 keV), of an amplitude simila r to that of the modu lation 
in the second channel (3-5 keV; cf. IClarkson et al. 20031 who 
report little or no variation below 3 keV). 

The most likely re ason for the discrep ancies between our 
analysis and that of IClarkson et al.l (120031) is the extended 
length of the lightcurve we were able to analy ze (about twice 
the length of that used in lOarkson et al. 2003) and our more 
precise definition of the edges of each individual period. The 
smoothing of super-orbital period length variation they ob- 
served is most likely at least partially an artifact of their use 
of a sliding Lomb-Scargle periodogram, and is similar to that 
we observed in Fig. [6] 

The question then arises as to why SMC X-l shows such 
extreme variations in its super-oribtal cycle length, while 
sources such as Her X-l and LMC X-4 do not. Several 
tentative suggestions have been made based upon radiatively 
driven warp models. In such models, the driving force of 
the warp is the flux from the central source, which scales 
as R~ 2 , being absorbed and reradiated from the outer disk, 
whose intercepting area scales as R 2 , thereby yielding a torque 
that increases linearly with distance. Disks can therefore 
become unstable to radiative driving of a warp if they are 
large enough (and the accretion efficiency is high enough; 
iPringld fl996l) . which is a function of the system's orbital 



period and secondary to primary m ass ratio. Numerical 
simulations dWiiers & Pringle 1 19991) and analytic estimates 
( Ogilvie & Dubusl l2001l) indicate that the Her X- 1 system is 
just large enough to have a stable warp. The analytic mod- 
els, however, also indicate that the, rather complex in terms 
of their eigenvalue behavior, underlying equations only ad- 
mit stable warp solutio ns for a finite range of disk sizes 
dOgilvie & Dubusll2001l) . Estimates are that SMC X-l may 
be near the upper limit of system sizes that are unstable to 
warping, and therefore may be more likely to exhibit chaotic 
warp modes. Our results are certainly consistent with that ex- 
pectation. 

Within the framework of a radiatively driven warp, the 
existence and stability of the warp modes are also depen- 
dent upon the radius at which matter is injected into the disk 
(IWijers & Pringlell999trOgilvie & Dubusl2001l) . Being more 
heavily wind-fed, this may be another way in which SMC X- 1 
differs from Her X- 1 . Perhaps variations in the wind speed of 
the SMC X-l secondary, leading to variations of the disk cir- 
cularization radius (equated with the mass injection radius in 
the majority of theoretical models) may account for the vary- 
ing cycle lengths. An ASM study like the one presented here, 
coupled with spectroscopic observations of the secondary's 
wind, may be an avenue for exploring such a possibility. 

The above speculations, however, do not account for 
the possibility of a fundamental "underlying clock" for the 
SMC X-l super-orb ital cycles. In their a nalysis of several X- 
ray binary systems iBovd & Sm ale (2004) find many systems 
for which the super-orbital excursion lengths are random inte- 
ger multiples of some fundamental period, (Cyg X-2, LMC X- 
3 and Cyg X-3). They also find, when plotting the lightcurves 
of the systems in X,X phase space, that each of these sources 
demonstrates circulation in the phase space about two rotation 
centers, creating a bimodal distribution in the length of super- 
o rbital variations. SMC X- l shows similar characteristics. 

Wiiers & Pringle (1999) speculated that for sufficiently 
large amplitude warps, the effective mass injection radius 
(i.e., where the incoming accretion flow interacts with and cir- 
cularizes into the inclined disk) could greatly fluctuate from 
one super-orbital cycle to the next, causing fluctuations in cy- 
cle duration that become tied to the orbital period. An avenue 
for exploring this possibility may be to model individual cy- 
cles of the ASM lightcurve, to see if one can discern a chang- 
ing warp amplitude with each cycle. The fact, however, that 
we did not find a dependence of cycle duration upon mean, 
minimum, or peak ASM flux, may argue against this possi- 
bility. 

To summarize, using the ASM lightcurve, we have de- 
fined a super-orbital ephemeris for SMC X-l for 1995-2006. 
It will now be possible to put pointed observations during 
this period in the conte xt of both the orbital ephemeris from 
Woidowski et al] d 1998b and the super-orbital ephemeris from 
this work. We have exploited the spectral information from 
the three channels of the ASM as much as possible and ex- 
tracted information from the spectrum in both the orbital and 
super-orbital profiles, to guide both further observations, as 
well as theories, e.g., of disk warping, which currently do not 
describe behavior as complex as that seen here in SMC X-l. 
It is also now necessary to perform multiple pointed observa- 
tions in order to observe a variety of super-orbital and orbital 
phase combinations and gather better statistics than provided 
by the data analyzed in this study. It might then be possible to 
analyze the super-orbital evolution of orbital profile features, 
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including the bounceback. With pointed observations we 
might also be able to better analyze the super-orbital ingress 
and egress and the spectral variation in the emission, to obtain 
more insight into the structure of the HMXB SMC X-l and 
other similar systems. 
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